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Abstract 

Diffraction, in the context of semiclassical mechanics, describes the man- 
ner in which quantum mechanics smooths over discontinuities in the classical 
mechanics. An important example is a billiard with sharp corners; its semi- 
classical quantisation requires the inclusion of diffractive periodic orbits in 
addition to classical periodic orbits. In this paper we construct the corre- 
sponding zeta function and apply it to a scattering problem which has only 
diffractive periodic orbits. We find that the resonances are accurately given 

by the zeros of the diffractive zeta function. 
PACS numbers: 0.320.+i, 03.65.Sq 
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Typeset using REVT^ 



Periodic orbit theory |l| encapsulates the duahty between local, classical information such 
as the periods, actions and stabilities of periodic orbits and global quantum information such 
as the density of states. Because of its reliance on classical mechanics, the theory encounters 
difficulties whenever the classical mechanics is singular. Examples of singularities include 
three-body collisions as in the Helium atom [Q; grazing conditions where some trajectories 
hit a smooth billiard surface while very close, parallel trajectories do not bouncing from 
a wedge where trajectories on one side of the vertex are reflected differently from those on 
the other side [^-^J; and, scattering from a point scatterer [^,|| or a magnetic flux line P] for 
which trajectories can not be continued through the discontinuities. In all of these examples, 
quantum mechanics smooths over the discontinuities through diffraction. 

In what follows, we study the third example - that of wedges. If one is interested in 
finding the trace of the Green function g{E) = TtG{E) (and hence the density of states 
through p{E) = —lmg{E)/n) of such a system, one must include the effect of not just 
classical periodic orbits but also so-called diffractive orbits These are paths which go 

directly into at least one vertex. Such paths obey classical mechanics everywhere but at the 
vertex. There one allows the path to enter and leave the vertex at any angle by assigning 
it an amplitude obtained by comparison with the exact solution of the quantum scattering 
problem [|l^ . 

This leads to the result that the contribution to the Green function of a diffractive path 
from point to via the vertex qy is approximately 

G(]^{qB, qA, k) ^ G^{qB, qy, k)d{e, e')G^{qv, q^, k). (1) 

Henceforth we assume a billiard system in two dimensions so that G'£(q2, qi, /c) = 
—iHQ^\k\q^2 — qi|)/4. The diffraction constant is 

,jn nf-^ ^ 4 sin (tt/ u) sm{9/u) sm{9'/u) 

^ ' ' V [cos (tt/z/) - cos {{e + e')lv)\ [cos (tt/z/) - cos {{e - e')lv)\ ■ ^ ' 

The angles Q and Q' are the incoming and outgoing angles relative to the same face of 
the wedge, although the choice of face is arbitrary. The wedge is parameterised by z/ = a/vr 



where a is the opening angle of the wedge. We have assumed Dirichlet boundary conditions 
on all surfaces. Note that d{6,6') = when a = Tc/n. For these special angles, we can 
continue any trajectory through the vertex by flipping the wedge n times to cover the plane 
0. In that event, the contribution to the Green function is not diffractive but geometric 
with a phase factor of (—1)" due to the n specular reflections all trajectories experience in 
the wedge. Eq. (|l]) is approximate and valid for the points A and B far from the vertex. The 
approximation breaks down when 6 ±6' = it, a condition which corresponds to the outgoing 
angle being directly on the border between a shadowed and an illuminated region as defined 
by the incoming angle. We return to this problem later. 

As with geometric orbits, the trace of G is a sum over periodic orbits. A stationary phase 
evaluation yields the contribution of a periodic diffractive orbit labelled 7 as [^,|| 

il i d ] 

^7(^) = -777 S n / ^f, } exp {i{kL.^ + n^vr - 3^^7r/4)}. (3) 
^fc [j=i ^StcH^j J 

The diffractive orbit has /x^ intersections with a vertex (diffractions) each with a correspond- 
ing diffraction constant d^j and reflections off straight hard walls. The total length of 
the orbit = ^-yj is the sum of the lengths of the diffractive legs along the orbit and 
is the length of the corresponding primitive orbit. If the orbit 7 is itself primitive then 
L-y = /-y. If 7 is the m'th repeat of some shorter orbit /3,then = Z^, = mZ/3, = mpp 
and fi^ = map. We then write 

9^ = 9p,m = -'^g (4) 

where 

= I n -Jm= \ exp {t{klf3 + p^n - Sa^n/A)}. (5) 
[j=i ^Snklpj J 

This follows because the contributions from the various diffractive legs in Eq. are mul- 
tiplicative so that Q/s^m can be factorised. This is not true for geometric (nondiffractive) 
periodic orbits. 



As with geometric orbits W^, we can organise the sum over diffractive orbits as a sum 



over the primitive diffractive orbits and a sum over the repetitions 



^d(^) =E £ = -4E^/3T^- (6) 

P m=l ^'^ /3 ^ 

We cast this as a logarithmic derivative by noting that ^ = ilj^tp — apti3/2k and recognising 
that the first term dominates in the semiclassical hmit. It follows that 

In addition to the diffractive orbits, one must also sum over the geometric (nondiffractive) 
periodic orbits so that in the logarithmic derivative we should also multiply by the contri- 
butions from the geometric orbits |^. In what follows, we assume that all periodic orbits 
are diffractive so that the poles of g{k) are the zeros of the zeta function |1T2|- [1^] 

C-\k)=Xl{l-t,). (8) 



For geometric orbits, this is evaluated using a cycle expansion [|T3|,|TJ]. Here the weights t/j 
are multiplicative so the zeta function is a finite polynomial conveniently represented as the 
determinant of a Markov graph |]TB|. 



It is instructive to consider a system which can be quantised solely in terms of periodic 
diffractive orbits, such as the geometry of Fig. 1. The classical mechanics consists of free 
motion followed by specular reflections off the sides of the wedges. The two vertices are 
sources of diffraction. The choice 71 = 72 has been used to study microwave waveguides 



and conduction in mesoscopic devices and is known to have at least one bound state ||T6|,|T7 . 
Unfortunately, in this case 9 + 9' = -k for the periodic orbits labelled B and B' in Fig. lb; 
Eq. (0) then diverges and the diffractive picture breaks down, as mentioned above. Instead, 
we consider 71 > 72. This is an open system with no bound states, only scattering resonances. 
As these are poles of g{k) and hence zeros of C^^(^) we can test the effectiveness of the theory 
in predicting them. 

In what follows, we consider the case 72 = for which the exact results are simple to 
obtain. (All the analytical results, however, are valid for 72 7^ 0). We define four cases: i) 
7i = 7r/2; ii) 71 > 7r/2; iii) 71 = 7r/2n (n > 2); and iv) all other values of 71. Cases i) and 



iii) differ from ii) and iv) respectively in that 71 corresponds to a special geometric angle 
and the lower vertex is not a source of diffraction. The reflection symmetry of the problem 
implies that all resonances are either even or odd. As a result, the zeta function factorises 
as (^^ = C+^C-^ and we determine and (Z^ separately. 

For cases i) and ii) there is only one periodic diffractive orbit which is labelled A in 
Fig. la. In the first case there is only one diffraction point while in the second case there 
are two. The weight of the periodic orbit in the two cases is found from Eq. (^) as tA = 



dAAGx.p{i(2kL + n / 4)} / Vl(hikL and tA = c/^ac^aa 6xp{2(2A;L + 7r/2)}/87rA;L respectively. 
The diffraction constants cIaa and d'^j^ refer to diffraction from the top and bottom vertices 
respectively and can be found from Eq. (0). The resonances are determined by = 1 which 
yields (to leading order) ||^ 
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i) kL^mT — 7i/8 In 



^IGn^rni — 7r/8) 



2 I dAA 

11) kL^n-n -'n/A--\xi[ — ^- — ' M . (9) 

2 \ «AA«AA / 

where n is a positive integer. Because the sole periodic orbit is on the symmetry axis, 



all of the resonances are predicted to be of even parity |18 . Unlike scattering systems 



quantised with geometric orbits p|,|5| jr3| , p!9[| , here there are no subleading resonances, a fact 



also observed in Refs. [^,0. This is a result of the multiplicativity of the Green function (p. 

For comparison, the exact resonances are found as follows. Defining polar coordinates 
with respect to the lower vertex, we expand the (unnormalisable) resonance wave function 



as 



ij{r,e) = ^anSm{an9)Jar,{kr) r<L 

n=l 

00 

= Y.br,sm{PJ)Hj;J{kr) r>L, (10) 

n=l 

where a„ = {2n + l)7r/27i and /5„ = mr/'-fi. This satisfies all the boundary conditions and 
we have restricted the discussion to even resonances. Demanding that the wavefunction and 
its normal derivative be continuous along the arc r = L gives the quantisation condition 
detM = where 



2 n 

M„™ = (-1)-+--^ wUkL) (11) 

and Wnm{z) is the Wronskian of Ja^{z) and Hp^{z). This relatively simple solution is a 
result of taking 72 = 0. 

The results for both the exact resonances and their semiclassical approximations are 
shown in Fig. 2a for 71 = 7r/2 and 71 = tt — 0.5. The weight Ia is greater by 0{\^kL) in the 
first case so that the resonances are not as unstable (i.e. not as far down in the complex k 
plane.) Although the agreement is very good, the theory is not very rich since it is based on 
just one periodic orbit. Of greater interest is a situation in which there are more diffractive 
orbits. As 71 is decreased, a pair of diffractive orbits is born each time 71 passes through a 
special angle 7r/2n. In particular, for '7r/4<7i<7r/2 there are three fundamental periodic 
orbits [Q, as sketched in Fig. lb, together with an infinite hierarchy of longer orbits labelled 
by their itineraries among the points A, B and B' . For the special case 71 = 7r/4, the lower 
vertex is not diffractive but rather induces two specular reflections. 

The weight of any orbit in which a letter is repeated can be expressed as the product 
of shorter weights (eg. tABB'B = tAsts's-) This means that the system can be written 
as a finite transfer matrix and the zeta function is the determinant of the corresponding 
Markov graph |20|. Due to the parity symmetry of the problem, we can consider just the 
right half of Fig. lb as the fundamental domain. Fig. 3 is the Markov graph of the system 
and shows all the ways of connecting points A and B. For example, the line marked as 
1 denotes starting at point B, going to the upper vertex and diffracting back. Symbols 
with a bar over them correspond to ending on the left half of the figure and then reflecting 
back onto the fundamental domain. For example, 2 denotes starting at B diffracting via 
the upper vertex to B' and then reflecting back onto B. Barred symbols contribute with a 
relative positive (negative) sign for the even (odd) resonances |T^. Symbol 5 is a so-called 
boundary orbit which lies on the border of the fundamental domain and contributes only 
to the even resonances. Each symbol carries its own weight and these combine to give the 
weights of the periodic orbits. To find an expression for the zeta function, we enumerate all 
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non- intersecting closed loops and non-intersecting products of closed loops to obtain |2G 



C±' = 1 - ti T ^2 - ^3(^4 ± h) - ^-^{h - hh T hh). (12) 

Since by symmetry = tg, we have 

= l-ti-t2-t5- 2^U + h{h + ^2) 
C' = l-(ti-t2). (13) 

We must still define the weights which appear in these formulas. Each symbol is com- 
posed of two segments leading from the corresponding nodes to the vertex connected by one 
diffraction at the vertex. We make use of this fact by separately defining quantities which 
contain the information about the segments and about the diffractions. The segment infor- 
mation is contained in ua and ub which are given by u\ = rf^^ exp{i{2kL — 3TT/A)}/87ikL for 



7i > Tc/A or u\ = exp{i2kL} /^/iGnkL for 71 = 7r/4, and u% = exp{i{2kH + n)} /^/iQnkH . 
These are simply the square roots of the weights of the fundamental periodic orbits shown 
in Fig. lb but without the phase and diffraction constant from the vertex. That information 
we quantify by defining the phase factor s = exp{— ■j37r/4} and four diffraction constants. 
There is a constant to diffract from A back to itself which for 72 = is d^A = 2. With 
similar notation we find (Iab = 2csc(7i/2 + 7r/4) and (Ibb = c^bb' = 1 + csc7i. (Note that 
by symmetry dij = dji.) The equality of dBB and dBB' is an artifact of choosing 72 = and 
is not a general result. We then have ti = sdBBU%, = sdBB'U%, = = ti = sdAsUA^B 
and ts = sdAA^A- 

We stress that the only difference between cases iii) and iv) is the form of ua and that 
the functional form ([T^) applies to both. However, as we will see, the numerical results are 
quite different. Due to the equahty oi dBB and dBB', the weights ti and are equal so 
there are no odd resonances. This can also be seen from the fact that Jp{z) and Hj)^\z) 
are independent functions so there can be no odd resonances which satisfy all the boundary 
conditions and match smoothly at r = L. 

We show the results for the even resonances for the cases 71 = 71/ 3 and 71 = tt/4 in 
Figs. 2b and 2c respectively. In the first case there is a diffraction at the lower vertex so 



that ua/ub = 0(1/ VkL). Therefore the dominant behaviour for is dictated by ti and 
^2 and the other terms represent small corrections. This is apparent in the figure where the 
gross behaviour is similar to Fig. 2a but with a decreasing, oscillatory correction coming 
from ua- In the second case ua is of the same order as ub so that ti, and ts are all of 
similar magnitude Then the resonances have no obvious pattern but are scattered around 
the complex k plane because of strong interferences among the various terms in (l^^. 

In both cases the semiclassical approximation is very accurate. In addition it gives us 
a better qualitative understanding of the spectrum than the exact numerical calculation 
based on Eq. ([To|) . Although we can not yet study the case 71 = 72 in detail, we can 



begin to understand why none of its bound states are odd [T^]. Recall that bound states 
correspond to real zeros of the zeta functions. However, the odd zeta function receives 
no contribution from the boundary orbit and also suffers from strong cancellation due to 
the relative negative sign between ti and t2, so only quite far from the real k axis are the 
magnitudes of the weights sufficiently large to allow for a zero. It is reasonable to suppose 
that this general condition continues to hold as 71 ^72. However, a complete semiclassical 
analysis of this case requires an understanding of the behaviour of the Green function on 
the border between shadowed and illuminated regions. 

Finally, we mention that the restriction to 71 > 7r/4 is not crucial. For smaller values of 
the angle, more diffractive orbits appear but the formalism above still applies and can be 
used to find the corresponding zeta function. 
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FIGURES 

FIG. 1. a) A configuration with only one diffr active periodic orbit, b) A configuration with an 
infinity of diffractive periodic orbits labelled by their itineraries among A, B and B' . 

FIG. 2. The resonances in the complex k plane, a) Results with only one diffractive periodic 
orbit. The upper set of points corresponds to 71 = 7r/2 and the lower set corresponds to 71 = tt— 0.5. 
For each set, the exact resonances are denoted with crosses and the semiclassical approximations 
with diamonds. With the same convention, b) shows the results for 71 = 7r/3 and c) shows the 
results for 71 = 7r/4. 



FIG. 3. The Markov graph used to compute the desymmetrised zeta functions. 
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